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TANGENT AND BERNOULLI NUMBERS RELATED TO MOTZKIN 
AND CATALAN NUMBERS BY MEANS OF NUMERICAL 

TRIANGLES 

JOSE LUIS ARREGUI 

Abstract. It is shown that Bernoulli numbers and tangent numbers (the derivatives 
of the tangent funcion at 0) can be obtained by means of easily defined triangles of 
numbers in several ways, some of them very similar to the Catalan triangle and a 
Motzkin-like triangle. Our starting point in order to show this is a new expression of 
C{2n) involving Motzkin paths. 



1. Introduction. Numerical triangles generated by matrices. 
Recall that the sequence of tangent numbers is defined as (tan(2"-i)(0)), with 

^ tan(2«-i)(0) 2„ , ^, , 

n=l ^ 

On the other hand, the sequence of Bernoulli numbers {B„) is defined by 

oo „ 

T^E^"^ (kl<27r). 
1 n\ 



z 



n=0 

The first values of Bernoulli numbers are 

Bf] = 1, Bi = — , Bo = — , Br = 0, Ba = , B^ = 0, Ba = — , 

' 2 6 30 42 

and in general B2n+i — for all n G N, while the signs of B{2n) alternate. 

From the relation between the exponential, sine and cosine functions, it results that 

4" (4" - 1) 



tan(2"-i)(0) = |S2„| 



2n 

Bernoulli numbers satisfy the recurrence relation 

Bn = r I , ] Bk 



k=0 



k 



These numbers are also related to Euler numbers {En), given by E2n+i = and 
E2n — (— where i?* are the secant numbers verifying 

secz=^ = y -^z^" (|z|<7r/2). 
cosz ^ (2n)! ^' ' ' ' 
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Actually, if we consider (6n) the sequence such that 
1 °° 

- (secz + tanz) = ^ (5„z"~^ (|z| < 7r/2) 

it turns out that 

(1) Vt-^ '—=Sn— (neN) 

^ ' (2fc + 1)« 2" ^ ' 

a result that goes back to Euler. For even values it follows the relation between Bernoulli 
numbers and the Riemann zeta function C, 

k=l ^ ' 

Standard references for all of this are Q and . 

A much more recent way to relate Euler and Bernoulli numbers is their combinatorial 
interpretation, due to R. C. Entringer, in terms of alternating permutations. See the last 
section for the details. 

In 1^ E. Calabi, F. Beukers and J. A. C. Kolk used a remarkable change of variables 
to reduce the series in (1) to the volume of a polytope in M", showing how it yields an 
elementary proof of (1). Very recently, N. D. Elkies has shown in ||^ how the same idea 
(which seems to be commonly attributed to Calabi) gives another proof of Entringer's 
theorem. In this paper we follow Calabi's argument in a different way, and although it 
applies only for even values of n (thus focusing on Bernoulli but not on Euler numbers) 
it allows us to find a nice connection between C,{1n) and Motzkin paths and numbers. In 
view of Entringer's result, a similar connection involves Catalan numbers as well. 

Let us first introduce one definition in order to provide an appropiate setting for our 
purpose: 

Let (A'^"^)„>i be a sequence of matrices (with numerical entries), each one with n 
rows and n + 1 columns. Let io,o = 1, = (ii,Oi^i,i) = ^^^^ and, for each n > 1, 

f — (f r. f , f n f \ — f 1 — 4(1) 4(2) . . . 4(n) 
We obtain a numerical triangle 



T — (^n.rn)o<7n<n — 



*0,0 














Hi 


^2,2 






^3,2 






^4,2 



^3,3 

^4,3 ti,- 



what we define as the triangle generated by (A(")). 

In particular, if A = {aij)i<ij is an infinite matrix and A(") is the submatrix formed 
by the first n rows and n + 1 columns, the triangle obtained in this way is the triangle 

generated by A. For instance, ^3 = (t3,07 ^3,1, ^3.2, ^3,3) = A^^"* A'-^'> A'-'^\ that is 




TANGENT, MOTZKIN AND CATALAN NUMBERS 3 

Example 1. If = 1 for all the the triangle generated hy A= (ay) is just 

tnm=n\ for all < TO < n. 



Example 2. Now let = < 



2 ifi=j=l, 
1 ifi=j-l, 
otherwise, 



and A = {ai 



(2 1 ■••\ 

1 •• 
1 •• 



7 



The triangle generated by A is given by tn,m = 2" that is 



1 

2 1 

4 2 1 

8 4 2 1 

16 8 4 2 1 



Example 3. For 



1 if i < i < i + 1, 
otherwise 



we have that 



{aij) 



generates the Pascal triangle 



/I 1 

110 

110 

1 1 

V 



10 10 



what means that tr. 



n 



TO 



for all n and to. Each number plus the one next to its left 



in the triangle gives the number below. 
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Example 4. Take the matrix in the previous example (say the Pascal matrix) and put 
1 in all the entries for the lower dominant as well, that is 

/I 1 •••\ 
1 



A 



Then the triangle generated by A begins 

1 

1 1 

2 2 
4 5 
9 12 
21 30 



1 

3 
9 

25 



1 

4 

14 



This is the so-called Motzkin triangle. Each number in the triangle is the sum of the one 
above and its (one or two) contiguous ones. 

The first column in the triangle gives the sequence of Motzkin numbers Mn- 

Example 5. If we fill with 1 all below the diagonal in the Pascal matrix we get 



A = 



/I 1 
1110 
11110 
11111 

V 



which generates 



1 

1 1 

2 2 
5 5 
14 14 



42 42 28 14 



that is known as Catalan triangle. Note that (if we set tn-\ = for all n) 

n 



k=m—l 



The first (second) column in the triangle gives the sequence of Catalan numbers Cn- It 
is well known that 

/2n — rn\ m + 1 
tn,m = I ) — — (0 < m < n) 



and thus C„ = 



2n\ 1 



n J n + 1 
(2n)! 



n J n + 1 n! (n + 1)! ' 



Both Motzkin and (much better known) Catalan numbers have a nice variety of re- 
alizations as the number of solutions to some combinatorial problem on n objects (see 
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|5j and [gj). In this article we will briefly explore one of them, obtaining a new one that 
links these numbers to permutation groups and then to Bernoulli numbers. 

For instance, we will obtain the following theorem: 
Theorem 1.1. Let A = (a^ ) the matrix given by 



-{ 




+ 1) 




-1 









ifi<3 


-1 


(2 


6 








•A 


2 


6 


12 







2 


6 


12 


20 




2 


6 


12 


20 30 




V- 








•7 



The first column of the triangle generated by A is then the sequence of tangent numbers 
(tan(2"-i)(G)) 

Note that A is the product of the matrix in example 5 and (on the right) the diagonal 
infinite matrix + l)(5y). 
The triangle begins 



1 








2 


6 






16 


48 


72 




272 


816 


1440 


1440 


7936 


23808 


44352 


57600 



where t„+i,m = m{m + 1) ^ tn,k- 

k—m — l 

We find it convenient to start with Motzkin paths and numbers, then showing their 
relation with Calabi's idea. 



2. Motzkin paths 
A Motzkin path of n-th order (or n steps) is a finite sequence 

n I 

A = (Ai, . . . , A„) e {-1, 0, 1}" such that ^ A^ = and ^ A^ > if / < n. 

The name Motzkin is due to the fact that the number of Motzkin paths of n steps is 
just the number M„ in example 4, so this is one of the realization of Motzkin numbers 
(surely the most common by now), and the term path comes after the usual visualitazion, 
in an obvious 1-1 correspondence, of any such sequence with a path that joins (0, 0) with 
(n, 0) in n steps -in the 'Lxl lattice-, each one by summing the vector (1, Aj), with the 
condition that no intermediate point lies below the horizontal axis. 

The empty sequence is considered as the only Motzkin path of steps, to give AIq = 1. 
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The only Motzkin path of 1 step is (0), and with 2 steps we have (0,0) and (1, —1). 
The four Motzkin paths of 3 steps are 

(0,0,0) /\ _ (1,-1,0) 

(0,1,-1) _ /\ / ^ \ (1,0,-1) 

and a 6 step Motzkin path is (1,0,1,— 1,0,— 1), visuaHzed as 

/ \ 

We wiU denote A4n the set of all Motzkin paths of n steps, being then M„ its cardi- 
nality. Ai — UJ^]^ Ain is the set of all Motzkin paths of any positive order. Being so 
easy to draw Motzkin paths, the reader can check that the first Motzkin numbers are 
1,1,2,4,9,21 (beginning by n = 0). 

The usual way to obtain Motzkin numbers is as in example 4. Another recursive 
formula is 

n-l 

M„+i = M„ + ^ AffeM„„fc„i, 

fe=0 

and in the next section we will recover the formula that best relates them to Catalan 
numbers. These and further results can be seen in ||l[ and ||^. 

To fix notations, let (A, p.) = (Ai, . . . , A„j , /ii, . . . , ^,12) if A = (Ai, . . . , A„j) and p. ~ 
(/ii, . . . , /ina)- Here p, could be the empty sequence, and (A, p) = A. 

The following proposition is then pictorially evident; it says that by removing flat 
steps, or by "flattening cusps" in a Motzkin path, we get another Motzkin path. 

Proposition 2.1. For any n G N 

(i) (X,0,p) G A4n+i if and only if (X,p) G A4„. 

(ii) (A, 1, -1, /i) G Mn+i if and only if (A, 0, p) G Mn- □ 
For each fc G N, let -Ife = (-1, . . . , -1), so (A, -1^) = (A, -1, . . . , -1). We allow 

k times k times 

A: = by defining — Iq as the empty sequence. In the same way, we define Ife and Ofc for 
any k G NU {0}. 

Now we define, for any p G {—1,0, l}*^, A4„;p = {(X,p) G Mn}, the set of Motzkin 
paths of n steps that "end in /2" , and M„.p is its cardinality. 

Note that Mn-j is empty, and then Mn = Mn;0 U Mn;-i- Moreover we have the 
following: 

Proposition 2.2. For every rt G N 

Mn = Mn-O U Mn:l-1 U Mn-Q,-1 U Mn-l,-l2 U Mn;Q,-l2 U Mn-1,-1^ U . . . 

(disjoint union), where the nonempty sets are exactly the first n sets in the list. 
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Proof. Note that Mn;-ik = Mn;i,-ik U Mn-fi,-ik U Mn;i-ik+i U . . . is empty if > 2n, 
since for any (A, —Ik) £ Mn wc have = Ai + • — h A„_fe — k < n — 2k. 

If n = 2k then (1^, — Ifc) € Mn-i-i^, while the same argument as before shows that 
Mn-fl-ik is empty. 

If n = 2fc + 1 then (1^,0, -Ife) G X„;o-ifc and (0, lfe,-lfe) £ Mna-i^- 
Finally, for n>2k (0„_2fc, Ifc, -Ifc) G and (0„_2fe-i, lfc,0, -Ife) G At„;o,-i;,- 

□ 



This decomposition of Mn allows us to prove the following "general" theorem, in the 
context of triangles generated by matrices. 

Theorem 2.3. Let f:M.^M.a function verifying the following recursion conditions 
(I) and (II), in terms of two matrices {bn,k)i<n,o<k o,nd {cn,k)i<n,k with bn,o 7^ for all 
n: 

(I) If (A, -Ifc) e Mn then /(A, 0, -U) = 5„,fe/(A, -U). 

(II) If (A, 0, -Ik-i) G Mn then /(A, 1^ -U) = c„,fc/(A, 0, -l^-i). 

Then the sequence (^n.o X^agai •^('^))n>i the first column of the numerical triangle 
f{0)T, where T is the triangle generated by the sequence of matrices (A'^"'), A^") being 
the n X (n + 1) matrix 

fbn.o c„,i 
bn,o bn,l 

bn^O bn,l Cn,2 
bn^O bn,l 6„,2 
bn,0 6„,i bn,2 C„,3 



V&n.O bn,l 6„,2 

Proof. Let (a;„,i,a;„,2, . • . ,Xn,n+i) be the vector given by 

( E /(^)' E /W' E /W' E •••) 



AeA1„+i;o 



AeA1„ + l;l,-l 



Aex„+i;o,-i 



Aex„+i;i,_i2 



I.e. 



E fW (0 < 2fc < n), 

AeA^n + l;0,-l;. 

a;„,2fe = E -^W (2<2fc<n+l). 

AeA^„+i;i,-ij. 

We claim that this vector equals f{0)A^^^A'-^^ ■ ■ ■ as we sec by induction on n: for 

n = 1 it is (/(O, 0), /(I, -1)) = (6i,o/(0), ci,i/(0)) = /(O)A(i). Assuming it true for n-1, 
note first that 

Xn,l= E /W= E /(A,0)= 5^ /(A,0) 

AeMn + l;0 (A,o)eA^n + i xem^, 

n 

= E ^",0/(^) = &n,0 E /(^) = ^n,oE^"-l'J' 

AeM„ AeAi„ j=i 
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For 1 < fc < n/2 we have 



AeA4„ + i;0,-ifc (A,-lfc)eM„ 

= bn,ki Y fw+ Y fw+ Y 

AeA^„;i,_i^. AGA^„;o,-i^. AeA^„;i,_i^,^j 

/ \ 





i>n,k 
\bn.k/ 



and finally, for 1 < A; < 



n + 1 



Ae7H„ + i;i,_i^ (A,0,-lfc_i)GA<„ 

^ c„,fe/(A,0,-lfc_i) =c„,fc ^ /(A) 

(A,0,-lfc_i)eA1„ AeA4„;o,-ij^_i 

/ \ 





V y 



so our claim is true. 



□ 



Our first application of this theorem is counting Motzkin paths. 

Corollary 2.4. The sequence oj Motzkin numbers (M„) is the first column of the triangle 
generated by A — (oij) with 

^ N 1 if i 'i'S odd and i > j — 1, 
I otherwise . 



Proof. Apply Theorem 2.3 to / = 1. 



□ 
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Note that 



A 



(I 


1 














1 





1 











1 





1 


1 








1 





1 





1 





1 





1 





1 


1 


1 





1 





1 






and the triangle (which is not the Motzkin triangle of example 4) begins 



1 
1 
2 
4 
9 
21 



1 
2 

5 

12 



1 

3 1 



Once the (n — l)-th row is given, we get the number tn^m as in Catalan triangle for odd 
m, but tn,„L = tn-i,m-i if "T- IS even. 

Following the proof of Theorem 2.3, note that the n-th row is just 

(Af„+1;0, Afn+l;l,-l, A^n + 1;0,-1, ^^n+l;l,-l2: ' ' ' ) ■ 

For instance, the third row is (4, 2, 2, 1), the cardinalities of 

A^4:o = {(0, 0, 0, 0), (1, -1, 0, 0), (1, 0, -1, 0), (0, 1, -1, 0)} , 

>(4;l,-l = {(0,0, 1,-1), (1,-1, 1,-1)}, 

A^4;0,-l = {(1,0, 0,-1), (0,1, 0,-1)}, 

A^4;l,-1. -{(1,1,-1,-1)}, 

which sum up 9 = M4 = M^.q = t^^. 



3. DYCK paths AND CATALAN NUMBERS. 

There is antoher easy decomposition of Motzkin paths that relates Motzkin and Catalan 
numbers, by counting how many Xj = there are in each A. For each n G N and 
A; = 0, 1, . . . , n, let 

Vn,k ^ {XeM; A has fc O's} 

and let Dn^k its cardinality. I?n,/c is then the number of Motzkin paths of n steps, k of 
which are horizontal. 

Since any Motzkin path has the same number of I's and — I's, the number of horizontal 
steps is of the same parity as n. Hence Dn.k = if n is even and k is odd or vice-versa. 

It is also obvious that — 1 for each n. 
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Vnfi is the set of the so-caUed Dyck paths, i.e. Motzkin paths with no horizontal steps. 
It is empty for any odd n, and then 

is the set of Dyck paths of any order. It is weU known that 

for aU n, where C„ is the Catalan number of example 5. See for instance [Q. We will 
find out a simple proof of this. 

With the same notations as in the previous section, we clearly have 

T^2nfi — Z^2n,0;l,-l U Z^2n,0;l,-l2 U ■ • • U 2?2ru0;l,-l„ 

(disjoint union). Denote the cardinality of each set by -D2n,0;i,-ifc > 0. 
Moreover, note that 

(A, -Ik) e I?2n,0:-lfc if and only if (A, 1, -Ifc+l) G I'2n+2,0:l,-U-+i , 

and 

T^2n,a:-lk = 2?2n,0;l,-lfc U T^2nfi\l,-lk + i U • ■ • U P2n,0;l,-1„ ■ 



Similarly to Theorem 2.3, we can formulate a "general" theorem: 



Theorem 3.1. Let f: Vq M a function verifying the following recursion condition, in 
terms of a matrix {an^k)i<n.o<k with a„.o =/= for all n: 
If (A, -Ife) e V2„fi then /(A, 1, -Ik+i) = an,fc/(A, -Ife)- 

Then the sequence (a„^o ■f(^))n>i the first column of the numerical triangle 

/(I,— 1)T, where T is the triangle generated by the sequence of matrices (A*^"^), A^"^ 
being the n x (n + 1) matrix 

( fln.o a„4 • • • \ 

a-nfi an,i an, 2 • • • 

an,l an,2 Ctn^S • • • 

O-n,0 On,! An, 2 On, 3 0,n,i ' ' ' 



\O'n,0 Cln,l On, 2 In, 3 On, 4 '•■ a-n,n) 



Proof. The proof is similar to that of Theorem |2.3| . We have to prove that 

Ae='C'2Ti + 2, 0:1,-1 -^£^2ti + 2,0;1,-12 ''^£'^2ti + 2,0; 1, - l^^i 

equals /(I, ••■ A'"), the first component being a„,o E /(^)- 

AeX'2„,o 



□ 



Take in particular f — 1. Then a„jt = 1 for all < ^ < and A is just the matrix of 
example 5. Theorem 3.1 says in this case that 



Cn = E] ^ ~ ^2ti,I 
A6l'2„,o 
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As in "our" Motzkin triangle, we also get that the numbers of the n-th row in the 
Catalan triangle are the number of Dyck paths of 2n steps finishing in — Ifc, for k = 
0, . . . , n. For instance, the third row is 

(5, 5, 3, 1) ~ {Dq q, Dq^q.^i, Dq q.^i^, Dq q.^i^) , 

the cardinalities of (listed backwards) 



^6,0;-l3 


- {(1,1,1,-1,-1, 










^6,0;-l2 


-P6A-13 U{(l,- 


1,1,1,-1, 


-1),(1,1, 


-1,1,-1, 


-1)}, 


2^6,0;-l 


-p6,0;-l2 U{(1,- 


1,1,-1,1, 


-1),(1,1, 


-1,-1,1, 


-1)}, 


2^6,0 


= ^6,0;-l • 











Our results concerning Bernoulli numbers are applications of Theorems |2.3| and BJJ 



Before seing that, we show another application of Theorem 2.3 
Theorem 3.2. Let 

fx 1/x ••■\ 



Aix) 



a; X 

X X 1/x 

X X X 

X X X 1/x 

X X X 



and let T{x) the triangle generated by A{x). Then its first column is given by the poly- 
nomials 



Dn,k 



Proof. Let z(X) the function that counts the zeroes in each A G A4, and let /(A) = x^^''^^. 
Note that 

/(A,l,-lfc) = i/(A,0,-lfc_i) and 
/(A,0,-lfc) =a;/(A,-lfc), 



so Teorem 2.3 applies, and the matrices A^") for each x are the corresponding submatrices 
of A{x). According to this, the first element of the n-th row is 



(A) 



fc=0 



□ 



Some remarks are in order: of course, if x = 1 it results the matrix of Corollary p.4| 
what is just saying that A/„ ~ J2k=Q ^n,k, obvious since {Vn,k)o<k<n give a decomposi- 
tion of A^„. 



On the other hand, note that from each Dyck path of P2m-2fe,o we get, by interlacing 
2k O's in all the (2") possible ways, this number of different paths in I?2m,2fc, and any 
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element in this set is obtained in this way from a unique path in I52m-2fe,o- It follows 
that 



Analogously 



_(2m + l\ _/2m+l 



These two facts give the most usual relation between Motzkin and Catalan numbers, 



n>0 



4. Calabi's argument and Motzkin paths. 

Let us briefly explain Calabi's argument reducing the calculation of C,{2n) to the volume 
of a polytope in R": 

For any natural number n > 2, by writing 



we get 



C(n)(l-2-")=f:^^. 



/c=0 



1 

Since — = / x^^dx, then 

2k + l Jo 



= ^ — T / (2^i^2 • • • Xnf''dm{x), 

^ - ^ t^r. J (0,1)" 



fe=0-^(0,l)'^ 

where of course x = {xi,. . . , a;„) and m stands for the Lebesgue measure. Then 

2'' 



C(n) 



2' 



T / y^XxiX2 ■ ■ ■ Xnf^dm{x) 

^ "'(0,1)" rrr; 



(0.^^- fc=0 

2" r 1 



2" — 1 J(o,l)" 1 ~ a^i^^i ' ■ ' 
Try now the change of variables 



dm{x) . 



, , /Sinui smu2 sinu„_i sinM^. 

{Xi,X2,...,Xn) = V{U) = ( , ,.. ., , ) . 

^008^2 cos Ms COSMn COS Ui ' 

sm u sill u 

Note that if < m, v < 7r/2 then — — ; > 0, and it is < 1 if and only if 

cost; sin(7r/2 — u) 

u < Tr/2 — V. Then, if we define 

A„ = {u> 0; Ml + U2 < 7i"/2, U2 + U3 < 7r/2, . . . , u„_i + < 7r/2, Un + Ui < 7r/2} 
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we have that <^(A„) C (0, 1)". It is not difficult to see that <^(A„) is actually (0, 1)" 
(moreover <^ is injective on A„), and hence this change of variables gives 



/ (0,1)" I — x\x2 - ■ -x"^ jj^^ 1 — tan^ Ui • • • tan^ u. 

Now, it is readily checked that 

det( Jv(u)) = 1 + (-1)"+^ tan^ m • • • tan^ u„, 
and thus, if n is even, it results 

/ 2~2 Trdmix) = [ dm = m(A„). 

i(o,i)- l~xixi---xl ^ ' Ja„ 

Take instead 

En = {Q < X & M"; Xi + X2 <l,X2+ Xi <1, . . . , Xn-l + X„ < l,Xn + Xi < 1} . 

Since m{An) = —m{En), we have shown that 

(2) C(2n) = ^^m(i?2„) 

for all n G N. Now let X2j = y2j and X2j-i = l—y2j-i for j = 1, . . . , n. This new change 
of variables gives 

where = {0 < y < 1; 2/2 < 2/i, 2/2 < 2/3, • • • , 2/2n < 2/2n-i, 2/2n < 2/i}- 

For n = 1, we have F2 = {{a,b) G M^; < 6 < a < 1}, so 771(^2) = 1/2 and 

C(2) = TrVe. 

Changing the coordinates order, we can re-write 

F2n = {{t,x) e (0,1)" X (0,1)"; h <Xi, ti<X2,...,tn < Xn, U < Xi}. 

Let ahb = min{a, 6} for all real numbers a and h. For each x G (0, 1)" the corresponding 
section of is 

F^ = {t; {t,x)eF2n} 

= (0, Xi A 0:2) X (0, a;2 A xs) X • • • X (0, A x„) x (0, a;„ A Xi), 



so 



m{F2n) = / m(Fs)dTO(x) = / {xi A X2) ■ ■ ■ {xn-i A Xn){xn A xi)dm{x), 
7(0,1)" -/(CI)" 

and we have obtained the following expression of ({^n) for each natural n: 



(3) C(2n) = / a^)dm{x) , 

* ~ J- ^(0,1)" 

where ^{x) = {xi A X2){x2 Axa)--- {xn-i A Xn){xn A xi). 



14 



J. L. ARREGUI 



For n = 2 we have ^(a, b) = {a A b)"^, and 

{a Abfdm{a,b) =2 {a A bfdm{a,b) = 2 b^dm{a,b) 

(0,1)2 J{0<6<o<l} J{0<b<a<l} 



= 2 f f b'^dbda =1 [ a^da = I 
Jo Jo 3 Jo 6 



so (3) says that C(4) = — • 



For n = 3 we can also write 

m{De) = 6 / ^{x)dm{x), 

J {0<X1<X2<X3<1} 

since any coordinate is compared to each other in {xi Ax2){x2 Ax3){x3 Axi). This makes 
it easy to see that m(i^6) = 1/15 and thus C(6) = 7r^/945. 

For n > 3 it does not suffice to consider the set {0 < xi < • • • < a;„ < 1}. 

Let Sn be the group of permutations of the set {1, 2, . . . , n}, for each n G N. 

Note that, in the integral in (3), we can ignore all points in (0,1)" with two equal 
coordinates, since 'm{H) = for any hypcrplanc H . Besides, if wc fix j and integrate 
in {x; Xj < Xi for all i ^ j, Xi ^ Xk for all i, k} then the result is independent from j, 
by simmetry. Choose j = n. For each x in the corresponding set there exists a unique 
permutation a S Sn-i such that 

<Xn< a;^(i) < a;<^(2) < • • • < a;<^(n-i) < 1- 
If Ecr is the set of such points for any given a, we have 



(4) / ^{x)dm{x) = n / ^{x)dm{x) = n / ^{x 



)dm{x) . 



Let E = Eii = {x; < Xn < xi < X2 < ■ ■ ■ < Xn-i < !}• 

The change of variables x i— > {x„(^i) , . . . , x„(^n-i) ) ^n) maps E onto E„-i , and then 



Ie 



^{x)dm{x) = / ^{x^(^i),...,x^(^n-i),Xn)dm{x) 
Je 

= / a;^(a=<7(i) A a;<,(2)) • • • (a;<,(„_2) A x„^ri-i))dm{x) . 
Je 



The last integral is easily expressed in terms of the ad hoc map 

U-2 5„^U~i{-l,0,ir 

defined as follows: for any n G N and a G S„+i, $(cr) is the family A = (Ai, . . . , A„) G 
{— 1, 0, 1}" such that, in the list 

0,cr(l),fT(2),...,cr(n),(T(n + l),0 

Xj + 1 is the number of neighbours of j that are greater than j. 

For example, if ct G S",.; is the given by 0, 2, 1, 4, 5, 3, then $(cr) = (1,-1,0,0): Ai + 1 = 
2 because 2 and 4 (the neighbours of 1) are both > 1; A2 + 1 =0 since 0, 1 < 2, and 
A3 + 1 = A4 + 1 = 1, for both 3 and 4 have 5 as a neighbour. 
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It is easy that — for slW A = ^{<^)- consider aU the couples of neighbours; 

appart from the first and the last (containing 0) each one adds +1 to ~ "-i 

there are n such couples. 

On the other hand it is clear that Ai > — 1 and A„ < 1, as in Motzkin paths. 

Definition 4.1. For each A G U,j5°^]^{ — 1, 0, 1}", J^(A) is the number of permutations a 
such that ^(ct) = A. 

Note that i^(A) > if and only if A is in the image of $. 

In the following Cj will denote (0, 0, . . . , 0, 1, 0, . . . ), and then Aicj = (Ai, A2, . . . , Aj ± 
1, . . . ). 

Lemma 4.2. For each A G { — 1,0, 1}" it holds: 

(i) 1^(1,0) ^2iy{X); 

(ii) v{\ -1) = 2 ^ v{\ - e,) + ^ z.(A - e,) . 

\j=0 \j=l 

Proof. Let a G Sn+i given by 0, cr(l), cr(2), . . . , cr(n), a{n + 1), 0. From a we obtain n + 2 
permutations r G Sn+2^ by interlacing n + 2 between two any neighbours. These r are 
different each other, and the (n + 2)! permutations of Sn+2 arise in this way by taking 
the (n + 1)! permutations of Sn+i- 

If r is so-derived from cr, with ^(cr) — A and $(t) = /2, then or Aj = /ij for j = 1, . . . , 71 
with = (this happens in two cases, setting n + 2 next to n + 1) or /ij = Aj + 1 for 
exactly one j, \j = /ij for the rest and jin+i = — 1- 

If Aj = — 1 then Hj = for each r that puts n + 2 next to j, and if Aj — then /ij = 1 
only if T sets n + 2 between j and its least neighbour. □ 



Theorem 4.3. The image of ^ is Ai, the set of Motzkin paths. 

Proof It is easy for n = I: Mi ^ {(0)}, ^2 = {(0, 1, 2, 0), (0, 2, 1, 0)} and i/(0) = 2. 
Assume that it is true for n. 



If A G and A = (/^,0), by Lemma 4.2 i^ip-) =/= 0, hence fi G Mn and then 



A G M.n+1 (by Lemma 2.1). If A = (/2, —1) Lemma says that, for some j, /2 — Cj has 
preimages in Sn+i, so /I — G A^„ and this easily implies that A G Mn+i (a simple 
draw may help). 

Conversely, if A G AAn+i and A = (/i,0) then /i G A^,,, and thus 2v{il) = i^(A) > 0, so 
A G $(5„+2). 

If A G Mn+1 with A = (/2, — 1) , by the definition of Motzkin paths we have j=i /"j = 1 
and Mi — fo'^ 1 < fc. Let I the maximum such that X]j=i P'j ~ ^- Sj<; /^j ~ 

it can't be ^ii — —1, so we can take — e;; Would it be /^j^i Mi = ^1 for some k > I, 
it should be X]j=fc+i A*j = 2, in contradiction with our election of L It follows that 



fj. — ei € Mn, and hence i'(/i — e;) > 0. Lemma 4.2 gives now that i^{fi, —1) = iy{X) > 



and hence A G <l'(5„+2). □ 
Let's get back to formula (4), where n > 3, and the integral 
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If A = ^'(ct) G Mn-2 then this integral is the same as 
_ ■ ■ ■ x^-^^L^"^^ dm{x) 

^n^2^^ ■ ■ ■ / x\^^^ / x\ dXn dxi ■ ■ ■ dXn-2 dXn-1 

Jo Jo 

111 1 



6n Ai + 5 Ai + A2 + 7 Ai H h A„_2 + (2n - 1) 



Definition 4.4. For eac/i A e {-1,0,1}", 

/9(A) = (Ai + 5)(Ai + A2 + 7) • ■ • (Ai + ■ • • + A„ + 2n + 3) . 

In this expression, each j-th factor {j > 1) is the previous one plus 3 (if A^ = 1), plus 
2 (if Aj = 0) or plus 1 (if A^ = -1). If A G AI then J^j = so the last factor is 2n + 3. 



Summarizing, we have that 

/ ^{x)dm{x) = ^ V —-^ 
Ao.D" 6 p($( 



^)) 



and then 



(5) / mdm{x)^l J2 ru^ 

which along with (2) gives, for every n > 3, 



Stop for a moment and compute C(6) and C(8): 
Ml ^ {(0)}, with j/(0) = 2 and p(0) = 5. Then 

TT^ 2 

C(6)- 



6(43 - 1) 5 945 ' 

M2 = {(0,0), (1,-1)}. From Lemma U i/(0, 0) = 2iy(0) = 4 and i^(l, -1) = i^(0) = 2. 
Besides p(0, 0) = 5 • 7 and p(l, —1) = 6-7, whence 

_8 /I r, _8 

C(8) 



6(44 - 1) ^5 • 7 6 • 7^ 9450 



Proposition 4.5. For any A £ U^]^{ — 1,0, 1}" iwe /lawe t/iat 
(i) v{X 0, -Ifc) = 2(fc + l)i/(A, -Ifc) /or each fc > 0, and 

fc + 1 - 

('iij v{\, 1, -Ife) = fc(fc + l)j^(A, -Ife-i) = — 2~^^^-' 0' ^Ife-i) for each k>l. 
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Proof. The proof follows by induction on fc, using Lemma 4.2. 

(i) If fc = the result is just (i) in Lemma 4.2. If fc > and we know that it is true for 
fc — 1, by (ii) in the lemma we can write 

z/(A,0,-lfe) = KA,0,-lfe_i,-l) 

= 2i^(A,-l,-lfc_i) + 2 ^ i^(A-e„0,-lfe_i)+ ^ i/(A - e, , 0, 



A, = l 



= 2i^(A, -Ife) + 2fc(2 J2 '^(^ - -lfe_i) +^2"^^' ^J' -Ifc-i)) 



Aj=0 



A, = l 



= 2i^(A,-lfc) + 2fci/(A,-lfc_i,-l) 
= 2(fc + l)KA,-lfe). 
(ii) If fc = 1 then 

:.(A,1,-1) = 2 ^ ^(A-e„l)+ ^ KA-e„l) + ^^(A,0) 



Ai=0 



Ai = l 



= zy(A,0) = 2i/(A) 

since 1) = for any /I. 

If fc > 1 and the result is true for fc — 1, then 

i^(A,l,-lfe) = i.(A,l,-lfc_i,-l) 

= 2 ^ v{X-ej,l,-lk-i)+ '^(A-ej,l,-lfe_i) + i^(A,0,-lfe_i) 

Aj=0 Aj=l 

= (fc - l)fc(2 ^(^ - e,. ~U-2) + ^ - ej, -lfe-2)) 

+ 2fcKA,-lfe-i) (by(i)) 
= (fc - l)fc z/(A, -lfe_2, -1) + 2fc //(A, -lfc_i) 
= (fc + l)fcz/(A,-lfe_i). 
The rest follows from (i): i^(A,0, — = 2fci/(A, — so 

z/(A,l,-lfc) = (fc + l)fc^i^(A,0,-lfc_i) = ^-iz/(A,0,-lfc_i). 

Corollary 4.6. Let A — (a.y) t/ie matrix given by 

2j if j is odd and i > j — 1 , 

J (j + 2) . . , . . -, 

' — if J IS even and i = ] — 1 , 

otherwise, 



I.e. 



A = 



/2 1 

2 4 

2 4 3/2 

2 4 6 

2 4 6 2 

V 



□ 
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Then, ifT— {tn,m) is the triangle generated by A, for all n 

tnfi = (ri + l)!. 



Proof. Proposition 4^ says that we can apply Theorem 2^ to f = i^, and conditions (I) 
and (II) there are satisfied with bn^k and Cn.k independent from n. It results that the 
corresponding A*^") matrices are the submatrices of A as stated, and since /(O) = 2 = 62.0 
we get 

which is just the cardinality of Sn+i, that is {n + 1)!. □ 



Proposition 4.7. (i) If k > and (A, — U-) G Mn then 

p{\, 0, -Ik) = p(X){2n + 5 - fc)(2n + 6 - fc) • • • (2?i + 4)(2n + 5), and 
p{X, -Ik) = p(A)(2n + 4 - k){2n + 5 - k) ■ ■ ■ {2n + 2){2n + 3) 
(just pCX) ifk = 0) 

(ii) If k> 1 and (A, 0, — 1^) G Mn, then 

p{X, 1, -Ik) = p(A)(2n + 5 - A:)(2n + 6-k)---{2n + 4){2n + 5). 

Proo/. If (A, 0, -Ife) e A^„+i then p(A, 0, -1^) = p{X)aia2 ■ ■ ■ ak+i, with ak+i = 2n + 5 
and aj+i — aj — 1 for each j, so 

p(A, 0, -Ik) = p(A)(2n + 5 - k){2n + 6 - k) ■ ■ ■ {2n + 4)(2n + 5). 

The rest is similar. □ 



Theorem 4.8. For each n >3, let 

(j + l)(4n-j + l) 



if j is odd and i < j — I , 



4)^^(,+2)(4n-,) = 
32 

^0 otherwise, 
and let A*") = (a|j'')i<,j<„_2, i<j<n-i- 

7^2" 4"-l 

C(2n) = - — bn , 

where 6„ is the first component of A^'^M^'^-' • • • A^"^ . 
Therefore 



B2n = (-1) 



n+l ^11 



2(4" - 1) 
tan(2"-i)(0) = bn. 
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^(4)^(^4 7/4 
U 7 



= 5 



^5 9/4 







.5 9 12/4 



/6 11/4 0\ 

6 11 

6 11 15/4 

\6 11 15/ 

/7 13/4 \ 

7 13 

7 13 18/4 

7 13 18 

V7 13 18 22/4/ 



and so on. 

Proof. Write formula (6) as 



_2n 

^(2") = 6(4^ _ 5: 



with /(A) = ^ for each XeM. 



From Propositions |4.5| and 4.7 it foUows that 
2(A: + l)(2n + 4-fc) 



/(A,0,-lfe) = 



(2n + 4)(2n + 5) 



/(A, -Ife) if (A, -Ifc) G Mn, and 



Hence / satisfies conditions / and // as in Theorem 2.3, with 

2(fc+l)(2n + 4-/s) 



(271 + 4j(2n + 5) 4 



2 

Note that /(O) = -, and then Theorem |2.3| gives, for any n > 3, 
5 



2n 



- J2 ,/(A) = (i(')i(")---i("',ei) 



where 4^'^+^) is the fc-th matrix in the statement of Theorem 2.3 and (•, •) is the usual 
scalar product. 

Let now 

^(„)^2n(2n+2)^(„).It 

is easily checked that A'-"-' are the matrices in the 

4 

statement, and we have now 

6-4'' 



so 



C(2n) = TT 



2n 4" ' 1 (^(3)^(4)... ^(n)^g^^ 



4" - 1 (2?i)! 



□ 
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5. Entringer's theorem and Catalan numbers. 



In order to prove Theorem 1.1 we will make use of Entringer's theorem about alternating 
(or zig-zag) permutations. 

A permutation a € Sn is said alternating if it is such that 

c^(j) < a-(j + 1) if and only if a{j + 1) > cr{j + 2) 

for j — 1, ... ,71 — 2, i.e. if a{j) is not a number between a{j — 1) and a{j + 1) for 
J = 2,...,n-1. 

Let a„ the number of alternating permutations. Let r £ Sn given by t(j) = n+l — j, 

i.e. 

T = n, n — l,n — 2,. ..,2,1. 
Then a t o a defines a bijection between alternating permutations a such that cr(l) < 
cr(2) and those such that cr(l) > cr(2), and thus the number of any of these is a„/2. 

Let then /3„ = Q!„/2, with Pq — /3i — 1. Entringer (see [Q) proved that these numbers 
give a combinatorial interpretation of tangent and secant numbers, namely 

secz = /3o + /32 y! + /34^ H and 

tanz = Piz + Ps^l + (3^^ + ■■■ 
3 5! 

for each z e C with \z\ < Tr/2. In particular 

/3(2n- 1) = tan(^"-i)(0) 

for aU n G N. 

Independently from Calabi's argument, in ||l2| R. Stanley obtained the tangent part by 
considering the polytopes in (1) and, as we mentioned in the introduction, N. D. Elkies 
(|^) has derived the result for both secant and tangent numbers starting from Calabi's 
idea. 

As for tangent numbers, there is the following relation with <I>: 

Proposition 5.1. Given n e N and a G S2n+i, ^(f) is a Dyek path if and only if a is 
alternating and (t(1) > cr(2). 

Proof. $((7) e 1^271,0 means that Xj = ±1 for all j, so in 

0, ct(1), cr(2), ct(3), . . . , a{2n), a{2n + 1), 

or two or none of the neighbours of each j are > j. This is exactly as saying that a is 
alternating and ct(1) > (j{2) (the fact that cr(2n) < <j{2n + 1) follows from them). □ 

Corollary 5.2. For each n e N 

tan(2"+i)(0) ^ KA). □ 

AGr'2,x,o 



Proof of Theorem 1.1 By Proposition |4.5| , for each fc > 

i^{X, 1, -Ife+i) = (fc + l)(fc + 2)z/(l, -Ik) 

whenever (A, — Ifc) € Vq. We can use Theorem |3.l| with f = v, and this together with 
Corollary ^^2] gives the result. □ 
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Let's see another application of Theorem 2.3, giving yet another way to obtain tangent 
(BernouUi) numbers by means of triangles. 

Theorem 5.3. Let T{x) he the triangle generated by the infinite matrix 



A{x) 



(2x l/x 

2a; 4a; 

2a; 4a; 3/2a; 

2a; 4a; 6a; 

2a; 4a; 6a; 2/a; 

V 



(so that A{1) is the matrix of Corollary ^.L ). Then the first element of the n-th row of 
T{x) is a polynomial Pn{x) of degree n such that Pn(l) = {n + 1)! CL'^id 

P„(0) =tan("+i)(0). 

Proof. Let f : M ^ M. given by fiX) — i'{\)x^^^\ with z(A) the number of O's in A, as 
in Theorem 3.2. By Proposition i.t we have 

/(A,0,-lfc) = 2(fc+l)a;/(A,-lfe), and 

/(A,l,-lfc) = ^/(A,0,-lfc_i), 
2a; 

and by Theorem |2.3| the first element in the n-th row of T{x) is 

n 

^ /(A) = KA)a;^(^^ = ^ "-^ =^ ^«(^) ' 
AeAi„ Aex„ fc=o 

where Un^k = ^(-^)- 

If n is odd then Pn {Q) = since 2?„^o is empty. 

If n is even, we have seen in Corollary O that P„(0) = tan("+i)(0). This holds for 
n = too, since Pq = I. 

Finally, = (n + 1)! as in Corollary O. □ 



Remark. Of course, the matrix ^(a;) is not defined for x — 0, but we certainly have 
that 

lim Pn{x) = tan("+i)(0). 

x—>-0 

If X is positive and small enough then tan("+^)(0) will equal the floor function applied 
to Pn{x). If we fix n and want to find tan'"+^^(0) it suffices to take x = l/{n + 1)!: we 
have seen that 

n 

P„(x)-tan("+i)(0) = 5]( KA))x^ 

k=l Ael?„,fc 

so if < a; < 1 and n is even 

n 

< Pn{x) - tan("+i)(0) <xY.i ''(^)) 

fe=l Aefn.fc 

<a;P„(l) = a;(n + l)!. 



and for x = 



1 

(n + 1)! 



we have < P„(x) - tan("+i)(0) < 1. 
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Entringer's paper |^] shows actually more than we have said so far: it provides a 
triangle to generate both tangent and secant numbers. Our last aim is to recall this, 
presenting it in the frame of our matrix-generated triangles idea. 

Let En^k denote the number of alternating permutations a G Sn+i such that cr(l) — 
k + 1 > (t(2), with Eo,o = 1. Note that Enfi = for n > 1, and i?„^„ is just the zig-zag 
number /3„. 

{En,k) are the so-called Entringcr numbers, and they form the Seidel-Entringer- Arnold 
triangle. The starting point for the proof of Entringer's theorem in [Q is the recurrence 

En+l,k+l = En+l,k + En^n-k [n > k > 0) , 

which gives itself a very simple algorithm to generate the triangle and thus the zig-zag 
numbers (see 

for 

more on this). This recurrence yields 

71 

(7) S„+i,fc+i = > ^ > 0) , 

j—n — k 

a fact that can be formulated as the following observation: if we remove the first (and 
trivial) column in the Seidel-Entringer- Arnold triangle, the resulting one is easily obtained 
by means of matrices: 

Proposition 5.4. For each n e N, let (^("^ — (a^ ) the n x {n+ 1) matrix given by 




if i + j l£ n , 

1 if i + j > n . 



Then the triangle (tn,m)o<m<n) generated by (A^"-*) is such that 

tn,m = £^n+i,m+i for all n, m . 
In particular tn.o is the zig-zag number /?„ for all n>0, and f„.„ = Pn+i- □ 



Note that A^"^ in the proposition is the submatrix of the first n rows and the last n-\-l 
columns of the infinite matrix 



and the triangle begins 



A = 



(■■■ 1 1^ 
••00111 
••01111 

) 



1 

1 1 

1 2 2 

2 4 5 5 

5 10 14 16 16 

16 32 46 56 61 61 
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Corollary 5.5. Let the infinite matrix 



A = 



/I 2 2 2 2 2 2 

1 2 3 3 3 3 3 

1 2 3 4 4 4 4 

1 2 3 4 5 5 5 

1 2 3 4 5 6 6 

1 2 3 4 5 6 7 

V 



and let A^"^ be the submatrix formed by the first n rows and the first n + 2 columns of 
A, for all n G N. Let then ti = (1) and 

in+l = t„+l,2, • • . , in+l,2n+l) = i^A^^""!) G N2"+^ (n G N) . 

Then, for every n = 1,2, . . . we have 

tn,i = p2n-2 (the u-th secant number), and 
tn,2n-i = p2n-i (the U-th tangent number). 

Proof. Just note that, if A^") are the matrices in the previous proposition, then 

^{2n-l) J{2n) ^ ^(2n-l) 

for each n G N, and thus the vectors t„ are the even rows in the triangle of the previous 
proposition. □ 



Proposition 5.6. Let A = (a^) the infinite matrix given by 



I.e. 









ifi 


is odd and 


i<3, 




ttij — 






if j 


is odd 


and 


j <i- 


f 


1, 








otherwise. 










n 


1 


1 


1 1 


1 


1 •• 


•\ 






1 





1 








• • 








1 





1 


1 1 


1 


1 •• 






A = 


1 





1 


1 





•• 








1 





1 


1 


1 


1 • • 








1 





1 


1 





1 ■• 

























Then the triangle {tn,m)n>m>o generated by A is such that 

tn,o = Pn+i (the zig-zag number) 
and each row tn is a permutation of the Entringer numbers {En+i^k)k=i,...,n+i- 
Proof. Using (7), it follows (by induction on n) that 

Tl 

tn,2k = -E„+i,„+i-fe if < fc < -, and 



tn,2k+l = if < < 



n-1 
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and then 



1 








1 


1 






2 


1 


2 




5 


2 


5 


4 


16 


5 


16 


10 


61 


16 


61 


32 



14 

56 46 
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